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Abstract 

Since  the  efficiency  of  fuel  cells  is  the  ratio  of  the  electrical  power  output  and  the  fuel  input,  it  is  a  function  of  power  density,  system  pressure, 
and  stoichiometric  ratios  of  hydrogen  and  oxygen.  Typically,  the  fuel  cell  efficiency  decreases  as  its  power  output  increases.  In  order  for  the  fuel 
cell  system  to  obtain  highly  efficient  operation  with  the  same  power  generation,  more  cells  and  other  auxiliaries  such  as  a  high-capacity  compressor 
system,  etc.  are  required.  In  other  words,  fuel  cell  efficiency  is  closely  related  to  fuel  cell  economics.  Therefore,  an  optimum  efficiency  should 
exist  and  should  result  in  the  definition  of  a  cost-effective  fuel  cell  system.  Using  a  multi-objective  optimization  technique,  the  sequential  quadratic 
programming  (SQP)  method,  the  efficiency  and  cost  of  a  fuel  cell  system  have  been  optimized  under  various  operating  conditions.  This  paper  has 
obtained  some  analytical  results  that  provide  a  useful  suggestion  for  the  design  of  a  cost-effective  fuel  cell  system  with  high  operation  efficiency. 
Published  by  Elsevier  B.V. 
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1.  Introduction 

A  fuel  cell  is  an  electro-chemical  energy  device  that  con¬ 
verts  the  chemical  energy  of  fuel  directly  into  electricity  and 
heat,  with  water  as  a  by-product  of  the  reaction.  As  a  renew¬ 
able  energy  source,  fuel  cells  are  widely  considered  one  of  the 
most  promising  energy  sources  because  of  their  high  energy 
efficiency,  extremely  low  emission  of  oxides  of  nitrogen  and 
sulfur,  and  very  low  noise,  as  well  as  the  cleanness  of  their 
energy  production.  Based  on  the  types  of  electrolytes  used,  they 
are  categorized  into  polymer  electrolyte  membrane  fuel  cells 
(PEMFCs),  solid  oxide  fuel  cells  (SOFCs),  phosphoric  acid  fuel 
cells  (PAFCs),  molten  carbonate  fuel  cells  (MCFCs),  and  direct 
methanol  fuel  cells  (DMFCs)  [1].  The  PEMFC  has  particularly 
attracted  more  attention  for  transportation  applications  because 
of  its  higher  power  density,  faster  start-up,  and  quick  response  to 
load  changes  than  other  fuel  cells  [2,3].  In  addition,  the  PEMFC 
operates  within  a  low  temperature  range  (60-80  °C)  and  has  a 
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relatively  simple  design  [1].  Due  to  these  multiple  advantages, 
the  PEMFC  has  become  the  best  candidate  for  an  alternative 
power  source  in  transportation  and  stationary  power  systems. 
To  commercialize  the  PEMFC,  cost  and  efficiency  need  to  be 
taken  into  account.  Thus,  achieving  an  optimal  PEMFC  system 
design  has  become  a  major  topic  in  recent  years.  For  stationary 
and  transportation  applications,  fuel  cells  are  required  to  achieve 
and  efficiency  equal  to  or  higher  than  40%  [3] .  However,  because 
the  efficiency  decreases  as  the  power  output  increases  [3],  as 
more  cells  request  power,  more  expenses  need  to  be  integrated 
to  achieve  a  high  efficiency  for  the  maximum  power  output.  To 
date,  although  many  techniques  of  optimization  of  fuel  cell  sys¬ 
tems  have  been  developed,  but  many  of  them  [5,6]  are  restricted 
to  only  one  optimization  objective,  such  as  performance  or  cost. 
However,  because  these  different  optimization  objectives  are 
coupled  or  affect  one  another,  considering  only  one  is  not  real¬ 
istic  or  practical.  Recently,  a  few  studies  of  multiple  objectives 
concerning  the  cost  and  performance  of  fuel  cell  systems  have 
been  reported  in  literature  [4] .  Xue  and  Dong  [4]  searched  for 
the  optimal  design  of  a  Ballard  fuel  cell  system  with  consider¬ 
ation  of  the  system’s  performance  and  cost.  In  their  work,  two 
system  parameters,  the  active  stack  intersection  area  and  the 
air  stoichiometric  ratio,  were  selected  as  the  design  variables 
in  the  joint  optimization.  Frangopoulos  and  Nakos  [7]  studied 
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Table  1 

Cell  voltage  parameters 


Parameter 

Value  and  definition 

N 

Cell  number 

E0 

Open-cell  voltage  (V) 

R 

Universal  gas  constant  (J  g-1  mol-1  K-1) 

T 

Temperature  of  the  fuel  cell  (K) 

F 

Faraday  constant  (C  mol-1) 

0 

Charge  transfer  coefficient 

i 

Output  current  density  (A  cm-2) 

io 

Exchange  current  density  (A  cm-2) 

in 

Internal  current  density  (A  cm-2) 

m  and  n 

Constants  in  the  mass  transfer  voltage 

r 

Area-specific  resistance  (kQ  cm-2) 

the  optimal  design  of  the  5-kW  PEM  fuel  cell  in  which  the  fuel 
cell  power  density  and  the  present  worth  of  the  life-cycle  cost 
of  the  system  were  used  as  the  objectives  of  optimization.  In 
our  study,  a  multi-objective  optimization  technique  is  applied  to 
optimize  two  objectives,  the  efficiency  and  the  cost  of  the  fuel 
cell  system  under  different  operating  conditions.  In  contrast  to 
previous  research  [4,7],  our  paper  defines  the  system  pressure, 
the  hydrogen  and  air  stoichiometric  ratios,  and  the  cell  voltage 
and  current  density  as  design  variables,  which  means  that  it  pro¬ 
vides  a  broader  insight  into  the  optimal  fuel  cell  design.  Under 
the  variations  in  these  variables,  the  structure  of  a  more  cost- 
effective  fuel  cell  system  with  a  high  efficiency  level  will  be 
determined  (Table  1). 

This  paper  is  organized  as  follows.  Section  2  gives  a  concept 
of  PEMFC  and  its  efficiency.  Section  3  addresses  the  design 
of  the  PEMFC  cost  model  for  the  optimization.  In  Section  4, 
the  multi-objective  optimization  for  the  PEMFC  is  presented. 
Section  5  provides  the  results  and  a  discussion  with  regard  to 
the  optimization.  Section  6  concludes  the  paper. 

2.  PEM  fuel  cell  and  efficiency 

A  PEMFC  consists  of  the  polymer  electrolyte  membrane 
sandwiched  between  two  electrodes  (anode  and  cathode).  In  the 
electrolyte,  only  ions  can  exit,  while  electrons  are  not  allowed  to 
go  through.  Thus,  the  flow  of  electrons  needs  a  path  like  an  exter¬ 
nal  circuit  from  the  anode  to  the  cathode  to  produce  electricity 
because  of  a  potential  difference  between  the  anode  and  cathode. 
The  overall  electro-chemical  reactions  for  a  PEM  fuel  cell  fed 
with  a  hydrogen-containing  anode  gas  and  an  oxygen-containing 
cathode  gas,  are  given  as  follows: 

Anode  :  2H2  4H+  +  4e_ 

Cathode  :  02  +  4H+  +  4e“  2H20 

Overall  :  2H2  +  02  2H20  +  electricity  +  heat 

In  practice,  a  5-kW  fuel  cell  stack,  such  as  a  Ballard  MK5-E 
PEMFC  stack,  uses  a  pressurized  hydrogen  tank  at  10  atm  and 
oxygen  taken  from  atmospheric  air  [10]. 

At  the  anode  side,  a  fuel  processor,  called  a  reformer,  which 
generates  hydrogen  through  reforming  methane  or  other  fuels 


Fig.  1.  Polarization  curve  (Ballard  Mark  V  PEMFC  at  70  °C)  [1]. 


like  natural  gas,  can  be  used  instead  of  the  pressurized  hydrogen 
tank.  A  pressure  regulator  and  purging  of  the  hydrogen  compo¬ 
nent  are  also  needed.  At  the  cathode  side,  there  is  an  air  supply 
system  that  contains  a  compressor,  air  filter,  and  air  flow  con¬ 
troller  to  maintain  the  oxygen  partial  pressure  [  1 , 1 2, 14] .  At  both 
sides,  a  humidifier  is  required  to  prevent  dehydration  of  the  fuel 
cell  membrane  [1,12,14].  In  addition,  a  heat  exchanger,  water 
tank,  water  separator,  and  pump  may  also  be  needed  for  water 
and  heat  management  in  the  FC  systems  [1,12,14]. 

To  produce  a  higher  voltage,  multiple  cells  must  be  connected 
in  series.  Typically,  a  single  cell  produces  a  voltage  between  0 
and  1  V  based  on  the  polarization  I—V  curve,  which  expresses 
the  relationship  between  stack  voltage  and  load  current  [1].  As 
shown  in  Fig.  1,  this  voltage  is  nonlinear  and  mainly  depends  on 
current  density,  cell  temperature,  reactant  partial  pressure,  and 
membrane  humidity  [1]  (Fig.  2). 

The  output  stack  voltage  Vst  [1]  is  defined  as  a  function  of  the 
stack  current,  reactant  partial  pressures,  fuel  cell  temperature, 
and  a  membrane  humidity,  as  follows: 

Vst  =  E  —  Vactivation  —  Vohmic  —  ^concentration  •  (1) 

In  above  equation, 


(Ballard  Mark  V  PEMFC  at  70  °C)  [1].) 

H2  Blower  and  Misc  (Values, 


4% 


Fig.  2.  Breakdown  in  stack  and  BOP  component  cost  contribution  for  an  80-kW 
direct  hydrogen  fuel  cell  system  [17]. 
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Table  2 

Ballard  Mark  V  PEMFC  coefficients  [14] 


Coefficients 

Values  (T in  °C) 

£oc  (V) 

1.05 

C(V) 

4.01  x  10"2  to  1.4  x  10“4r 

r  (k£2  cm-2) 

4.77  x  10"4  to  3.32  x  10"6r 

m(V);T>  39°C 

1.1  x  10“4  to  1.2  x  io~6r 

m  (V);  T<  39°C 

3.3  x  10“3  to  8.2  x  10“5r 

n  (cm2  mA-1) 

8.0  x  10“3 

E  =  N[E0  +  (RT/2F)  ln((/fc2  /7\7)/l%Oc)]  is  the  ther- 
modynamic  potential  of  the  cell  or  reversible  voltage  based  on 
the  Nemst  equation  [1].  Vactivation  is  the  voltage  loss  due  to  the 
rate  of  reactions  on  the  surface  of  the  electrodes.  V0hmic  is  the 
ohmic  voltage  drop  from  the  resistances  of  proton  flow  in  the 
electrolyte.  VCOncentration  is  the  voltage  loss  from  the  reduction 
in  concentration  gases  or  the  transport  of  a  mass  of  oxygen  and 
hydrogen.  Their  equations  are  given  as  follows: 


_  RT 

^activation  =  N-  —  In 

2a  F 


=  NCn 


(2) 


Vohm  =  Nir  (3) 

^concentration  =  Nm  exp(m).  (4) 

In  Eq.  (1),  Pu 2,  Po2,  andJPH2oc  are  the  partial  pressures  of 
hydrogen,  oxygen,  and  water.  Subscript  4c’  means  that  the  water 
partial  pressure  is  vented  from  the  cathode  side. 

The  concentration  voltage  loss  is  simplified  by  neglecting  in, 
which  is  a  very  small  value  compared  with  i,  and  combining  i0 
of  the  Tafel  equation,  with  NA  ln(/)  [1].  In  practice,  the  Ballard 
mark  V  fuel  cell  voltage  per  cell  is  described  in  Eq.  (5),  which 
has  the  specific  coefficients  given  by  Table  2. 


Vc  =  E0 c  —  ri  —  C  ln(/)  —  m  exp (ni). 


(5) 


The  detail  explanation  of  each  voltage  loss  can  be  found  in 
Ref.  [1].  In  the  Nernst  equation,  the  ideal  standard  potential  E0 
for  a  PEMFC  is  1.229  V  with  a  liquid  water  product,  or  1.18  V 
with  a  gaseous  water  product  [8].  Under  the  assumption  that 
pressure  on  both  the  cathode  and  the  anode  is  approximately  the 
same,  the  Nernst  equation  is  transferred  into  a  function  of  the 
system  pressure  Rsys  [1,9]  given  as  follows. 


RT 

E  =  E0-\ - In 

0  2  F 


RT 

+  lnCAys) 


(6) 


where  a ,  and  8  are  constants  depending  on  the  molar  masses 
and  concentrations  of  H2,  O2,  and  H2O.  Each  partial  pressure 
can  be  expressed  by  these  constants  and  the  system  pressures. 


Ph2  —  Oi  Psys 

EO2  =  fiEsys  •  (2) 

^H20  =  SPsys 


Assuming  that  a,  /3 ,  and  8  are  constants,  Eq.  (9)  shows  that 
the  EMF  of  a  fuel  cell  is  increased  due  to  the  system  pressure 
Ps ys.  Hence,  Psys  can  be  used  as  one  of  the  optimization  design 
variables.  For  the  multi-objective  optimization,  the  specification 
of  the  fuel  cell  stacks  has  to  be  identified  in  advance,  and  then 
each  optimization  model  can  be  derived. 

First,  the  fuel  cell  efficiency  optimization  model  is  derived 
based  on  Ref.  [10],  and  the  output  power  of  the  fuel  cell  system 
is  described  by  the  following  equations. 

Pfcs  =  2*stack  —  Pprs-  (8) 

P stack  =  NiVcA  =  50  kW.  (9) 


Pprs  —  +  R0th- 


^comp  —  Op 


OmOmt 


Psys  A 

~P^J 


0.286 


—  1  I  m. 


m  =  3.57  x  10  7  x  A. air  x  i  x  A  x  N kgs  1 , 


(10) 

(ID 

(12) 


where  Rfcs  is  the  net  power  of  the  fuel  cell  system  and  Pstack  is 
the  stack  output  power.  The  parasitic  power  consumed  by  the 
compressor  is  Rcomp  and  the  others  is  Roth-  Even  though  Roth 
was  assumed  to  be  a  constant  of  2kW  in  Ref.  [10]  based  on 
62.5  kW-rated  stack  power,  here  it  is  5%  of  the  nominal  power 
out,  50  kW,  due  to  the  inclusion  of  unexpected  power  consump¬ 
tion,.  Thus,  R0th  is  assumed  to  be  2.5  kW.  The  flow  rate  of  air 
m  is  related  to  the  air  stoichimetry,  the  cell  current  density,  and 
the  active  cell  area.  Before  proceeding  to  build  the  efficiency 
optimization  model,  let  us  consider  how  to  determine  the  opti¬ 
mal  cell  number  and  cell  area.  Once  optimal  current  density  and 
cell  voltage  have  been  decided,  the  total  active  cell  area  ( N  x  A) 
can  be  calculated  by  using  optimal  power  density,  the  product 
of  Vc  and  /,  i.e.,  we  can  decide  on  a  number  of  cells,  as  long  as  a 
single  active  cell  area  is  given.  The  system  pressure  of  the  fuel 
cell  is  always  higher  than  the  atmospheric  pressure  in  a  certain 
range  because  the  compressor  cannot  provide  a  pressure  under 
the  atmospheric  pressure.  According  to  Ref.  [11],  Rsys  must  be 
0.02  MPa  higher  than  the  inlet  air  pressure  Rin,  so  it  is  included 
as  one  of  constraints  in  the  optimization.  Thus,  if  using  the  lower 
heating  value  (LHV),  the  fuel  cell  efficiency  optimization  model 
is  obtained  to  achieve  the  maximum  efficiency  of  the  fuel  cell 
inEq.  (13). 


max  ?7fc(Rsys>  ^H2>  ^ air?  Vc,  /)  — 


Vc^f(Rstack  Pprs) 
1.25Rstack 


S.t. 


Rsys  >  0.12  MPa 
-  1 

2- air  >  1  -  (13) 

Vc  >  ov 
i  >  0  mA  cm-2 


where  Uf  is  the  fuel  utilization  rate,  which  is  the  reverse  of  the 
hydrogen  stoichiometric  ratio  [12].  The  air  stoichiometric  ratio 
must  be  over  the  minimum  limit  in  order  to  prevent  the  depletion 
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of  oxygen  at  this  minimum  limit.  The  hydrogen  stoichiomet¬ 
ric  ratio  also  is  greater  than  1  unless  it  runs  in  the  hydrogen 
dead-ended  mode  [12].  Normally,  a  higher  air  and  hydrogen  sto¬ 
ichiometric  ratio  is  preferred  in  low  power  ranges.  The  ranges 
of  cell  voltage  and  current  density  will  be  based  on  the  V-I 
polarization  curve.  Since  the  cell  voltage  is  a  function  of  the  cell 
current  density,  as  seen  in  Eq.  (5),  we  can  select  four  optimization 
parameters:  system  pressure,  air  and  hydrogen  stoichiometric 
ratios,  and  cell  current  density.  By  using  four  optimization 
parameters,  the  optimization  model  has  been  built.  In  the  fol¬ 
lowing  section,  the  fuel  cell  cost  optimization  model  will  be 
presented. 

3.  PEM  fuel  cell  cost  model 

For  this  analysis,  we  are  particularly  interested  in  the  small 
and  middle-sized  fuel  cell  systems.  We  use  a  50-kW  PEM  fuel 
cell  system  for  transportation  as  the  example  of  our  study.  For 
the  fuel  cell  cost  model,  the  cost  of  fuel  cell  stack  and  balance  of 
plant  (BOP)  components  for  water,  thermal,  and  fuel  manage¬ 
ment  were  assessed.  Due  to  lack  of  the  latest  data  about  hydrogen 
storage,  power  electronics,  electric  drive  motors,  hybrid  batteries 
for  PEM  fuel  cell  systems,  the  fuel  storage  and  fuel  genera¬ 
tion  components  were  excluded  in  this  study.  The  target  cost 
is  the  cost  of  the  fuel  cell  stack  and  the  BOP  system,  given  as 
follows: 


Table  4 

Fuel  cell  stack  cost  parameters 


Production  volume 

Cost  parameter, 

A  (US$  m-2) 

Cost  parameter, 

B  (US$) 

100 

811.77 

1311.30 

1,000 

722.54 

363.33 

10,000 

454.45 

428.51 

30,000 

329.24 

405.79 

60,000 

312.26 

160.98 

The  parameter  A  is  the  power-dependent  term  with  regard  to 
US$  m-2  of  the  membrane  area,  and  the  parameter  B  is  the  fixed 
cost  for  the  fuel  cell  stack. 

In  Ref.  [14],  the  annual  fuel  cell  degradation  is  assumed 
to  be  a  6%  year-1  drop,  and  the  planned  fuel  cell  lifetime  is 
assumed  to  be  87,600  h,  10  years.  The  platinum  loading,  Lp, 
is  defined  as  0.4  mg  cm-2,  and  the  cost  of  platinum,  Cp,  is  US$ 
1160  troy-1  oz-1.  Although  this  Cst i  (US $)  model  is  used  in  Ref. 
[7],  it  is  not  logically  understandable  because  many  constants 
and  parameters  (M,  A,  and  B)  are  involved  in  Eq.  (15)  without 
justification.  Thus,  a  more  reasonable  stack  cost  model  Cst2  is 
chosen  from  Refs.  [17,18]  in  our  study.  The  Cst2  (US$  kW-1)  is 
described  as  follows: 

r  -  -  Cm  +  Ce  +  Cb  +  CPt  +  C°  i 


C  =  Cst  +  Crop, 


(14) 


Cpt  —  Cwp t  x  Tpt 


(17) 


where  Cs t  is  the  cost  of  the  fuel  cell  stack  and  Crop  is  the  cost 
of  the  balance  of  plants.  For  Cst,  currently  two  types  of  fuel 
cell  stack  cost  models  are  available  in  Refs.  [13,16].  One  is 
represented  by  the  following  equation  [13]. 


Csti  =  M 


105.4  17.56LpCp 


10 


380 


PG(\+df 

Pd 


+  B 


(15) 


where  M  is  a  fixed-cost  markup  (1.1  default);  A  and  B  the  cost 
parameter  that  depends  on  production;  volume  (see  Table  3);  Lp 
the  fuel  cell  platinum  loading  for  both  electrodes  (mg  cm-2);  Cp 
the  cost  of  platinum  (US$troy-1  oz-1);  Pq  the  fuel  cell  gross 
dc  peak  power  (kW);  P&  the  fuel  cell  power  density  (W  cm-2);  d 
the  annual  fuel  cell  degradation  (%  year-1);  and  A  is  the  planned 
fuel  cell  lifetime  (years). 


Table  3 

Specification  of  the  fuel  cell  system  based  on  Ref.  [1] 


P  =  10  X  Vc  X  i.  (18) 

where  Cst2  is  the  fuel  cell  stack  cost  per  kW  (US$kW-1);  Cm 
the  membrane  cost  (US$  m-2);  Ce  the  electrode  cost  (US$  m-2); 
Cb  the  bi-polar  plates  cost  (US$  m-2);  Cpt  the  cost  of  platinum 
catalyst  loading  (US$m-2);  Cwpt  the  weight  of  platinum  cata¬ 
lyst  loading  (g  m  -2x  Fpt  the  unit  cost  of  platinum  (US$g 
C0  is  the  cost  of  peripheral  materials  (US$m-2)  that  include 
end  plates,  plastic  frame,  and  thrust  volts;  Ca  the  assembly  cost 
(US$  kW-1);  Vc  the  cell  voltage;  and  i  is  the  cell  current  density 
(A  cm-2)  (Table  4). 

In  Refs.  [17,18]  the  cost  of  each  component  of  a  50-kW  PEM 
fuel  cell  stack  was  estimated  based  on  an  automatic  production 
line  with  an  annual  production  capacity  of  18,000  vehicles. 

Table  5  shows  the  specific  costs  of  components  in  the  PEM 
fuel  cell  stack. 

As  for  Crop*  the  cost  model  can  be  found  in  Ref.  [13]  - 
Crop,  including  the  air  blower,  humidification,  radiator,  stain¬ 
less  pump,  iron  pump,  control  electronics,  actuation,  piping,  and 


Items 

Specification 

Nominal  power  output  (kW) 

50 

Stack  temperature  (K) 

353  (80 °C) 

Inlet  H2/air  humidity  (%) 

100 

Cell  open  voltage,  E0  (V) 

1.05 

Entry  air  temperature,  Tc  (K) 

288  (15  °C) 

Specific  heat  constant,  cp  (JK-1  kg-1) 

1004 

Compressor  efficiency,  rjc 

0.75 

Compressor  connecting  efficiency,  r]m 

0.85 

Inlet  pressure,  P\n  (Pa) 

105 

Table  5 

Specific  costs  for  components  in  PEM  fuel  cell  stack  [17,18] 


Components 

Cost 

Nafion  membrane  (US$  m-2) 

500 

Platinum;  2-4  gm-2  (US$m-2) 

32-64 

Electrode;  max.  0.8  mm  for  single  cell  (US$  m-2) 

177 

Bi-polar  plate;  max.  4  mm  (US$  m-2) 

1650 

Peripheral  parts  (US$  m-2) 

15.6 

Assembly  (US$kW_1) 

7.7 
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valves  -  which  is  approximated  by  a  quadratic  equation  in  the 
fuel  cell  output  power  and  varies  with  the  production  volume, 
as  shown  the  following  equations  [13]. 


For  100  production  units: 

Ca  =  3343.5  +  39.942/>g  -  0.0454Pg. 

(17a) 

For  10,000  production  units: 

Ca  =  2980.2  +  35.654Pg  -  0.0422 P5. 

(17b) 

Unfortunately,  these  Eqs.  (17a,b)  are  valid  for  stationary  PEM 
fuel  cells  and  were  published  in  1999.  Instead  of  using  these  two 
equations,  the  most  updated  data  in  Ref.  [17]  are  used  in  our 
study,  in  which  Crop  is  estimated  to  be  34%  of  fuel  cell  system 
cost,  as  Cst,  including  assembly,  is  assumed  to  be  contributing 
approximately  66%  of  fuel  cell  system  cost.  Even  though  the 
breakdown  of  the  fuel  cell  system  is  for  an  80-kW  direct  hydro¬ 
gen  system,  the  same  breakdown  is  possible  to  use  for  a  50-kW 
fuel  cell  system  because  the  cost  analysis  is  only  integrated  in 
the  stack  and  BOP  costs. 

So,  in  order  to  build  the  cost  model  sharing  the  same  opti¬ 
mization  parameters,  the  cost  model  is  likely  to  be  a  function 
of  the  efficiency,  as  in  Eq.  (21),  which  is  able  to  investigate  the 
impact  of  the  cell  voltage,  Vc,  and  current  density,  i,  as  well  as 
other  optimization  parameters  on  the  fuel  cell  cost.  In  our  study, 
the  costs  of  the  fuel  cell  stack  and  BOP  are  considered,  and  the 
maximum  fuel  cells  system  cost  is  obtained  as  the  maximum  effi¬ 
ciency  is  achieved.  Thus,  this  cost  optimization  model  can  share 
the  same  optimization  parameters  with  the  system  efficiency 
optimization  model,  and  each  feasible  range  of  the  parameters 
will  be  used  as  a  constraint  of  this  optimization  problem.  The 
base  production  volume  is  chosen  as  18,000  units,  which  can  be 
a  mass  production  for  the  mobile  50-kW  PEM  fuel  cell  system. 

min  CpcC^sys,  7,h2  5  ^air?  Vc,  0  =  (Cst  +  Crop)  x  rjfc 
s.t. 

Psys  >  0.12  MPa 

7-h2  >  1 

7- air  >  1  •  (21) 

Vc  >  o  V 
i  >  0  mA  cm-2 

As  explained  in  Section  2,  the  design  variables  of  the  cost 
model  could  be  reduced  to  the  cell  current  density. 

In  the  following  section,  the  multi-objective  optimization 
will  be  presented  with  consideration  of  both  efficiency  and  cost 
optimization. 

4.  Multi-objective  optimization  for  PEMFC 

According  to  the  above  considerations,  the  multi- objective 
optimization  problem  is  formulated  as  follows: 


,u  ,  ,  .x  /  n  VcWf(Pstack  —  Pprs) 

min  ^fcC^sys?  ^H2  5  7. air,  l)  —  (—  1) 

I  stack 

minCpc^sys,  7,h2  »  ^air>  0  =  (Cst  +  Crop)  X  7?fc 
s.t. 

Psys  >  0.12  MPa 

^H2  >  1 

.  (21) 

7- air  >  1 

i  >  0  mA  cm-2 

In  the  multi-objective  optimization  problem,  both  objective 
functions  have  to  be  minimized  simultaneously.  The  objectives 
in  such  a  problem  are  often  in  conflict  with  each  other.  From 
Eq.  (21),  when  the  efficiency  is  increased,  the  cost  increases 
as  well.  Thus,  during  the  optimization  process  for  such  a  prob¬ 
lem,  there  is  no  single  optimum  solution  to  improve  all  of  the 
objectives,  which  means  that  a  number  of  solutions  that  are  all 
optimal  can  exist.  As  the  solution  to  the  multi-objective  problem 
is  a  set  of  points  that  represent  the  best  trade-offs  between  the 
objective  functions  for  each  solution,  there  is  no  way  to  further 
improve  an  objective  function  without  worsening  at  least  one 
other  one.  Such  points  are  called  Pareto-optimal  points  or  non¬ 
inferior  points.  The  set  of  all  the  Pareto-optimal  points  is  called 
the  Pareto-optimal  set  or  the  Pareto  frontier.  Here,  the  MATLAB 
optimization  toolbox  [16]  for  multi-objective  optimization  prob¬ 
lems  is  used  to  find  the  Pareto  frontier  solution  set.  MATLAB  has 
two  functions  to  solve  multi-objective  problems:  fminimax  and 
fgoalattain.  Even  though  both  methods  use  the  popular  nonlinear 
programming  algorithm,  a  sequential  quadratic  programming 
(SQP)  fminimax  method  is  more  appropriate  for  our  optimiza¬ 
tion  than  the  fgoalattain  method  because  the  latter  method  is 
more  complicated  than  the  fminimax  due  to  the  weighting  coef¬ 
ficients.  The  general  form  of  fminimax  method  is: 

minmax{/i,  fi, . . . ,  fm} 

x  f 

such  that 
Ax  <b 

^eq-^  =  bQq 

C(x)  <  0  ,  (22) 

Ceq(x)  =  0 
Up  ^  x  <  Up 

where  v  is  the  design  variable  vector,  f\  ,/2, . .  .,/m  are  the  objec¬ 
tive  functions,  matrix  A  and  vector  b  are  the  coefficients  of  the 
linear  inequality  constraints,  matrix  Aeq  and  vector  beq  are  the 
coefficients  of  the  linear  equality  constraints,  C  contains  the  non¬ 
linear  inequality  constraints,  Ceq  contains  the  nonlinear  equality 
constraints,  and  Lp  and  Up  are  the  lower  and  upper  bounds, 
respectively.  In  order  to  search  for  the  optimal  design  value  x,  the 
fminimax  method  iteratively  minimizes  the  worst-case  value  (or 
maximum)  of  the  objective  functions  subject  to  the  constraints. 
The  advantage  of  this  method  is  to  easily  find  the  optimum  design 
from  an  arbitrary  initial  design  point.  Furthermore,  less  func¬ 
tion  and  gradient  evaluation  are  required  compared  with  other 
constrained  nonlinear  optimization  methods.  However,  the  main 
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disadvantages  of  both  methods,  fminimax  and  fgoalattain,  are 
that  the  objective  functions  must  be  continuous  and  that  each 
method  has  a  limitation  in  searching  for  global  solutions. 

5.  Results  and  discussion 

In  this  section,  the  fminimax  method  is  executed  to  solve 
the  multi-objective  optimization.  To  avoid  the  unrealistic  design 
criterion,  the  upper  bounds  of  the  system  pressure  and  the  air 
and  hydrogen  stoichiometric  ratios  are  specified  at  10,  10  and 
10  MPa,  respectively.  According  to  the  polarization  I-V  curve 
in  Fig.  1,  the  cell  current  density  will  lie  within  the  range  of 
0-1  A  cm-2,  which  is  used  as  one  of  bound  limits  for  the  opti¬ 
mization.  The  bound  limits  of  the  design  variables  are  given  as 
follows: 

0.12  MPa  <  Psys  <  10  MPa 
1  <  A,h9  <  10 

“  .  (23) 

1  <  Aair  <10 

0  <  /  <1  A  cm-2 

With  various  initial  conditions  of  the  design  parameters,  Ps ys, 
Ah2,  ^air,  and  i,  corresponding  trade-offs  (Pareto)  solutions  are 
obtained.  So  to  speak,  these  solutions  are  in  the  Pareto  set;  that 
is,  as  one  objective  is  improved  in  the  set  the  other  is  worsened. 

For  simplicity,  the  initial  conditions  can  be  described  in  the 
column  vector,  such  as  [PSys>  ^h2>  ^air,  *]•  If  an  arbitrary  ini¬ 
tial  condition  is  defined  as  the  vector,  In  1  =  [0.12  MPa,  2,  2, 
800  mA  cm-2],  Fig.  3  shows  the  trade-off  solution  based  on  the 
vector  In  1 .  In  Fig.  3,  as  the  efficiency  is  improved  in  the  set,  the 
cost  is  increased  as  well.  In  changing  the  initial  condition,  this 
Pareto  frontier  will  be  changed  because  the  fminimax  method 
will  find  a  local  solution  in  the  changed  initial  condition.  First, 
when  In  1  is  changed  to  In  2  =  [0.24  MPa,  2,  2,  800  mA  cm-2]  - 
that  is,  the  system  pressure  becomes  two  times  larger  than  the 
one  for  In  1  -  Fig.  4  shows  how  this  change  affects  the  opti¬ 
mization  of  the  fuel  cell.  As  seen  in  Fig.  4,  the  higher  system 
pressure  is  able  to  achieve  a  cost-effective  and  high-performance 
model  compared  with  the  model  given  by  In  1 .  For  instance, 


Efficiency 

Fig.  3.  Pareto  frontier  based  on  the  In  1.  In  1  =  [0.12  MPa,  2,  2,  800mA  cm-2]. 
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Fig.  4.  Pareto  frontier  change  from  In  1  to  In  2.  In  2  =  [0.24  MPa,  2,  2, 
800  mA  cm-2]. 


Trade-off  (Pareto)  Solution 


Efficiency 


in  Fig.  4,  at  the  efficiency  0.45,  the  In  2  Pareto  frontier  cor¬ 
responds  to  about  420US$kW-1,  whereas  in  In  1  the  Pareto 
frontier  is  almost  700US$kW-1  than  mean  Ini  is  less  eco¬ 
nomic  condition  than  the  In  2.  However,  in  the  case  of  much 
higher  pressure,  the  change  from  In  1  to  In  3  =  [0.36  MPa,  2, 
2,  800  mA  cm-2],  as  shown  in  Fig.  5,  although  a  better  opti¬ 
mum  model  than  In  2  is  obtained,  the  efficiency  range  of  In  3 
becomes  unrealistic  because  the  fuel  cell  system  efficiency  is 
not  normally  greater  than  0.6  [1].  Hence,  the  initial  condition 
In  2  is  more  recommendable  than  In  3.  Second,  if  the  initial 
hydrogen  stoichiomertic  ratio  is  changed  from  2  to  1.5,  such  as 
In 4  =  [0.12  MPa,  1.5,  1,  800  mA  cm-2],  the  less  cost-effective 
and  poorer  performance  model  is  found  in  Fig.  6.  Hence,  around 
the  hydrogenstoichiomertic  ratio,  2  is  more  preferable  than  1.5. 

In  the  case  that  the  air  stoichiometric  ratio  is  changed  from  2 
to  1.5  as  the  In5  =  [0.12MPa,  1.5,  1,  800mA cm-2],  the  more 
cost-  and  efficiency-effective  model  is  achieved  in  Fig.  7.  How¬ 
ever,  if  it  keeps  decreasing  to  1 ,  the  efficiency  is  not  applicable 
to  the  real  system  as  seen  in  Fig.  8.  Hence,  the  ratio  around  1.5 
is  more  preferable  than  2  and  1.  For  the  current  density,  as  it 
decreases,  the  cost-  and  efficiency-effective  model  is  achieved. 
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Fig.  5.  Pareto  frontier  change  from  Ini  to  In 3.  In 3  =  [0.36 MPa,  2,  2, 
800mA/cm2]. 
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Fig.  6.  Pareto  frontier  change  from  Ini  to  In4.  In 4  =  [0.12 MPa,  1.5,  2, 
800  mA  cm-2]. 
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Fig.  7.  Pareto  frontier  change  from  Ini  to  In5.  In5  =  [0.12MPa,  2,  1.5, 
800  mA  cm-2]. 

Through  trial  and  error,  the  optimal  current  density  is  approxi¬ 
mately  estimated  to  be  450  mA  cm-2.  Comparing  this  with  In  1, 
In  7  provides  a  better  performance-  and  cost-effective  model, 
as  shown  in  Fig.  9.  With  the  recommended  current  density  of 
450  mA  cm-2  and  the  V-I  polarization  curve  shown  in  Fig.  1, 


Trade-off  (Pareto)Solution 


Trade-off  (Pareto)  Solution 


Fig.  9.  Pareto  frontier  change  from  Ini  to  In 7.  In 7  =  [0.12 MPa,  2,  2, 
450  mA  cm-2]. 


the  optimal  cell  voltage  can  be  calculated  to  be  0.72  V,  and  the 
power  density  will  be  3.2  kW  m-2.  If  a  50-kW  rate  power  out¬ 
put  is  selected,  then  the  total  active  cell  area  is  15.625  m2,  which 
means  the  stack  will  need  to  contain  174  layers  of  a  single  cell 
with  30  cm  x  30  cm  active  cell  area.  Thus,  the  optimal  current 
density  allows  us  to  determine  the  total  active  cell  area  and  even 
to  provide  the  information  about  the  cell  number  required  for 
the  target  stack  if  the  a  single  cell  area  is  given. 

Figs.  4-11  indicates  that  the  change  of  each  design  variable 
has  a  severe  impact  on  the  cost  and  efficiency  of  the  fuel  cell. 
Because  the  current  density,  in  particular,  is  closely  associated 
with  the  fuel  cell  area,  it  has  a  more  direct  effect  on  the  fuel  cell 
cost  and  efficiency  than  any  of  the  other  variables. 

When  selecting  the  preferable  initial  condition, 
In 8  =  [0.24 MPa,  2,  1.5,  450mA cm-2],  based  on  the  above 
discussion,  the  Pareto  solution  is  achieved  in  Fig.  10,  but  it  is 
definitely  not  applicable  to  the  real  system  due  to  an  impracti- 
cally  high  efficiency.  Therefore,  In  8  must  be  adjusted  such  that 
the  Parteto  solution  lies  in  a  realistic  range.  As  In  8  =  [0.24  MPa, 
2,  1.5,  450mA cm-2]  is  adjusted  to  In 9  =  [0.24 MPa,  1.75, 
1.5,  450mA cm-2],  this  condition  can  be  used  for  the  design 


Trade-off  (Pareto)Solution 


Fig.  8.  Pareto  frontier  change  from  Ini  to  In 6.  In 6  =  [0.12 MPa,  2,  1, 
800  mA  cm-2]. 
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Fig.  10.  Pareto  frontier  change  from  Ini  to  In 8.  In 8  =  [0.24 MPa,  2,  1.5, 
450  mA  cm-2]. 
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Efficiency 

Fig.  11.  Pareto  frontier  change  from  In  8  to  In  9.  In  9  =  [0.24  MPa,  1.75,  1.5, 
450  mA  cm-2]. 

of  an  economic  and  high-performance  fuel  cell.  Any  initial 
condition  can  be  chosen  as  long  as  it  is  within  the  bound 
limits  and  the  corresponding  Pareto  solutions  are  applicable  in 
practice. 

6.  Conclusion 

In  the  paper,  a  joint  optimization  model  of  fuel  cell  system 
efficiency  and  cost  is  proposed.  A  multi-objective  optimization 
technique,  the  sequential  quadratic  programming  method,  has 
been  applied  to  investigate  the  impact  of  the  variations  of  ini¬ 
tial  conditions  on  the  efficiency  and  cost  of  a  fuel  cell  system. 
Although  the  study  shows  that  the  change  of  current  density  is 
more  closely  related  to  the  fuel  cell  efficiency  and  cost  than  to 
any  other  variables,  the  system  pressure,  the  hydrogen  and  air 


stoichiometric  ratios,  and  the  current  density  must  be  appropri¬ 
ately  selected  for  the  optimal  design  because  they  also  largely 
affect  the  fuel  cell  system  and  cost,  as  well,  in  the  our  study.  Our 
work  presents  a  way  to  determine  the  optimal  design  regarding 
the  fuel  cell  efficiency  and  cost  aspects  simultaneously. 
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